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Abstract: Drought is one of the most common natural events that have a great negative impact on agriculture and water 
resources. Recently, a new index for drought assessment and monitoring is presented called Reconnaissance Drought Index 
(RDI). RDI is calculated based on precipitation and potential evapotranspiration. In this study, indices ofSPI (Standardized 
Precipitation Index) and RDI were calculated using 30 years meteorological data (1982 to 2012) for Kohgilouye and 
Boyerahmad Province located in iran. The negative values show the dry years of the historical record. It can be seen that 
the most sever drought occurred during the year 2007- 2010. Correlation analysis is performed to identify differences of the 
SPI and RDI in different time scales. Results show the presented correlation coefficients increase in the longer time scales 
such as (9 and 12 months). The Map of Spatial distribution of drought severity for 2010 obtained based on the SPI and RDI 
by means of geo statistics methods in the study area. It can be observed that East part of the watershed lie in the moderately 
class and the rest of area is generally characterized by severe drought. In this paper however it was shown that RDI can be 
used for monitoring purposes and to a certain extent for short period drought forecasting. 
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I. Introduction 

Drought as an environmental disaster is associated with a deficit of water resources over a large geographical area, 
which extends for a significant period of time (Rossi, 2000). It occurs in areas with high and low rainfall and all climate 
conditions. The different drought indices are presented for drought monitoring. They are useful tools for decision makers to 
identify weather abnormal conditions for a region (Wilhite et al., 2000). 

There are many literatures on the drought evaluation by using different indices, models and water balance 
simulations (Jain et al., 2010). There are different indices for drought monitoring, such as: Palmer Drought Severity Index 
(PDSI), Crop Moisture Index (CMI), Surface Water Supply Index (SWSI), Percent of Normal Index (PN), Standardized 
Precipitation Index (SPI) (Mishra and Singh, 2010). 

Along the various indices for meteorological drought monitoring, SPI were widely accepted and used (Shamsnia 
and Pirmoradian, 2009). The SPI was presented by McKee and his colleagues at Colorado State University . SPI is 
completely related to probability. This Index use precipitation as the most effective parameter in the calculation of the 
drought severity. The SPI is normally distributed so it can be applied to monitor wet as well as dry periods. During the last 
decade the SPI has become very popular due to its low data requirements. This index can be useful for drought monitoring 
and forecasting (Shamsnia et al., 2009). Recently, a new index for drought assessment and monitoring is presented called 
Reconnaissance Drought Index (Tsakiris et al., 2007). RDI is calculated based on precipitation and potential 
evapotranspiration. Precipitation alone cannot show the impact of drought on agricultural production and vegetation. 
Applying both of the P and PET in drought severity calculation and monitoring, increases the validity of the results (Tsakiris 
et al., 2007). Also, this index has a strong correlation with SPI (Asadi Zarch et al., 2011). The RDI was used in Greece 
(Tigkas, 2008), Cyprus (Pashiardis and Michaelides, 2008), Malta (Borg, 2009) and Iran (Khalili et al., 2011) for monitoring 
or analysis of historical droughts. To achieve this objective, we used Geo statistical tools to mapping. Geostatistical Analyst 
creates a continuous surface by using sample points taken at different locations. Using interpolation techniques are 
frequently used in recent years to mapping drought risk. 

The Study Area 

Drought is a common natural disaster in many countries as well as Iran. Previous studies show that Iran is exposing 
drought with different severities. To evaluate drought, Kohgilouye and Boyerahmad Province in the southwest of Iran were 
selected. This province is located in the southwest part of Iran, at 49° 57' to 51° 42 E longitude and 30° 9' to 31° 32' N 
latitude, with an area of 26416 Km2. Elevation in the study area ranges from 500 m to 4400 m. The climate is characterized 
by dry and warm summers, and cool wet winters. The annual mean of precipitation for the province ranges from 400 to 
1300mm. Most of the annual average precipitation falls as rain during winter. The Location of study areas and monitoring 
station are shown in Figs 1 . 

Data used 

Rainfall data were collected from 17 meteorological stations complied by the Meteorology Organization of Iran and 
Kohgilouye and Boyerahmad Regional Water Authority. After, annual rainfalls data during 30 years from 1982-2012 were 
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chosen to later analysis. Missing data was computed by use of relationship between elevation and rainfall. To this order, 
linear regression was calculated in Spss software. 

II. Method 

Calculation of meteorological indices 

The Standardized Precipitation Index 

The standardized Precipitation Index (SPI) is one of the indices that were presented for drought monitoring. It is 
calculated in short-term (3, 6 and 9 months) and long-term (12, 24 and 48 months) periods. In any time scale, the SPI mean 
may reach zero in a location and its variance become equal to 1. Using the SPI, quantitative definition of drought can be 
established for each time scale. A drought event for time scale i is defined here as a period in which the SPI is continuously 
negative and the SPI reaches a value of -1.0 or less. The drought begins when the SPI first falls below zero and ends with the 
positive value (Mckee et al., 1993). The SPI is calculated by using the following equation: 

P. — P 

SPI = (1) 

Where P t and P are the precipitation and average precipitation of selected period, S is standard deviation of precipitation. 
The Reconnaissance Drought Index (RDI) 

The Reconnaissance Drought Index (RDI) is calculated in three stages: Initial value of RDI (a 0 ), normalized RDI (RDIn) and 
standardized RDI (RDIst). Initial value may be calculated for each month, seasons (3 -month, 4-month, etc.) or hydrological 
year. The a 0 is calculated by using the following equation (Tsakiris et al., 2007; Asadi Zarch et al., 2011) 

yl2 p 

Where Pz/and PETy are the precipitation and potential evapotranspiration of the jth month of the ith hydrological year. 
Hydrological year is starting from October in Iran. N is the total number of years of the available data. 

A second step, the Normalized RDI(RDI n )is computed using the following equation for each year, in which it is evident that 
the parameter a 0 is the arithmetic mean of Rvalues (Tsakiris et al., 2007; Asadi Zarch et al., 2011). 

RDjf^-X (3) 
a 

The third step, the Standardized RDI(RDI S ), is computed following a similar procedure to the one that is used for the 
calculation of the SPI: The equation for the Standardized RDI is: 

(o _ yi ~ yk 



Where y k is the ln(a 0 ; y^) is its arithmetic mean of y k , and d y k is its standard deviation. 

The RDI is based on the ratio of two aggregated quantities which are precipitation and potential evapotranspiration. It can be 
estimated for all time scales. However 3, 6, 9 and 12 month are suggested since they are more useful for comparing different 
locations (Tsakiris and Vangelis, 2005). 

Estimation of potential evapotranspiration for drought monitoring using RDI 

Several methods have been suggested for PET estimation. The Penman-Monteith (PM) equation is one of the 
standard methods of calculating PET by using meteorological data. It is remarkable that the international scientific 
community has accepted this equation as the exact method (Jabloun and Sahli, 2008). Although PM equation in different 
climates has shown positive results, but because of the full meteorological data requirement such as minimum and maximum 
air temperature, minimum and maximum relative humidity, solar radiation and wind speed, its widespread use is limited. In 
some studies for drought monitoring using RDI, alternative methods such as Thornthwaite method and Blaney-Criddle 
equation has been used for PET estimation (Kanellou et al., 2010). Therefore in this research was used the Calibrated 
Thornthwaite method (Ahmadi and Fooladmand, 2008) for Kohgilouye and Boyerahmad Province. The Calibrated 
Thornthwaite equation is : 

T 

PET = 16(10 -^-) a (5) 

Where PET is mm month" 1 , Tejf is the effective monthly temperature (°C). / is a thermal index imposed by the local normal 
climatic temperature regime, and the exponent a is a function of /. 
Effective temperature is calculated as (Camargo et al., 1999): 

Teff = K{T ave +A) (6) 

^ 1 max T m i n (7) 
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Where Tmin, Tmaxand Teffaxt the minimum, maximum and effective monthly temperature (°C) respectively. Kis the calibrated 
coefficient. 

Time Scales 

In this study, Monthly SPI and RDI were calculated using the monthly precipitation and evapotranspiration data of 
17 stations in Kohgilouye and Boyerahmad Province. SPI and RDI values were calculated in the period of hydrological years 
from 1982 to 2012 for the time scales of 1, 3, 6, 9 and 12 months. Because of data limitations, SPI with time scales longer 
than 24 months may be unreliable. The RDI st behaves similar to the SPI. Drought classification of these indices is shown in 
Table 1. The values from -0.5 to -0.99 and +0.5 to +0.99 considered as nominal class of mild drought and slightly wet 
(Shamsnia and Pirmoradian, 2009). 

III. Results and discussions 

Correlation analysis between RDI and SPI 

The correlation coefficients (r) between SPI and RDI for each area and each time scales are described in Table 2. 
Results show the presented correlation coefficients increase in the longer time scales such as (9 and 12 months). The 
maximum value of r between SPI and RDI was obtained for 12 -month time scale and the minimum value for 1 -month time 
scale in all areas (Table 2). In time Scale 9 and 12 months in all stations, correlation coefficient (r) is more than 0.9. Because 
SPI index is related to precipitation parameters for drought assessment, the results show that precipitation parameters are 
more effective in drought with a longer time scale. By contrast, the effect of potential evapotraspiration parameters that were 
used in the RDI index will be reduced by increasing the time scale. The wet, drought and normal periods have more 
fluctuation on the short time scales such as 1, 3 and 6 months, because of the low number of effective months and effects of 
other meteorological parameters on the severity of drought. In other words, changes of dry and wet periods on short time 
scales, in addition to precipitation, depend on the evapotraspiration and the other weather parameters. 

Identification of Drought Events 

Fig 2 shows the annual SPI and RDI in the yasuj station for the period 1982 to 2012. The negative values show the 
dry years of the historical record. It can be seen that the most sever drought occurred during the year 2007- 2010. 

Drought Maps 

Drought maps were produced using the estimated annual values of SPI and RDI. Also, By means of kriging method 
spatial severity of drought was mapping. Fig 3 and 4 shows Spatial distribution map of SPI and RDI in 2010 in the study 
area. It can be observed that north and west parts of the area lie in the extremely class and the rest of area is generally 
characterized by severely and moderately drought. Such spatial drought map can be useful to forecast drought in the further 
studies. In most cases, the results from both indices are similar. However, in some cases, depending on the local conditions, 
the results from RDI are more sensitive. 

IV. Conclusions 

SPI and RDI were calculated in a period of 30 years data (1982 to 2012) for Kohgilouye and Boyerahmad Province. 
Correlation analysis is performed to identify the differences of the standardized precipitation index (SPI) and the 
reconnaissance drought index (RDI). From the presentation of some popular indices for drought assessment it seems that SPI 
is becoming the most widely used index. This is probably due to its simplicity, universality and to its least data demanding 
nature. In this paper however it was shown that water deficit cannot be estimated only on the input (e.g. precipitation) but 
also on the output variable (water consumption). Based on the logic a new index, RDI was purposed using data of two 
determinants, precipitation and potential evapotranspiration. It can be more effectively associated with hydrological and 
agricultural drought. Also it is an ideal index to study the effects of climate instability conditions. Due to its easy calculation 
RDI can be used for monitoring purposes and to a certain extent for short period drought forecasting. 
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Table 1 . Drought classification according to the SPI and RDIst values 



SPI and RDLt value 


Category 


2 or more 


Extremely wet 


+1.5 to +1.99 


Severely wet 


+1 to +1.49 


Moderately wet 


+0.5 to +0.99 


Slightly wet 


-0.49 to +0.49 


Normal 


-0.5 to -0.99 


Mild drought 


-1 to -1.49 


Moderately drought 


-1.5 to -1.99 


severely drought 


-2 and less 


Extremely drought 



Table 2. Correlation coefficient (r) of the SPI and RDI in different time scales 



Row 


Station 


Time Scales (Month) 






1 


3 


6 


9 


12 


1 


Yasuj 


0.838 


0.965 


0.979 


0.990 


0.992 


2 


Dashtroom 


0.521 


0.746 


0.899 


0.948 


0.985 


3 


Dehnoo 


0.640 


0.831 


0.954 


0.977 


0.991 


4 


Patave 


0.532 


0.764 


0.842 


0.971 


0.989 


5 


Abchirak 


0.615 


0.856 


0.941 


0.983 


0.997 


6 


Nazmakan 


0.654 


0.832 


0.932 


0.955 


0.984 


7 


Tolchega 


0.543 


0.683 


0.742 


0.886 


0.979 


8 


Ghaleraesi 


0.725 


0.846 


0.921 


0.964 


0.993 


9 


Bibihakime 


0.626 


0.754 


0.878 


0.943 


0.996 


10 


Likak 


0.537 


0.622 


0.774 


0.865 


0.947 




Legend 

-t climatology station "^^h==^^b=^^^ 

| | province boundary 1 : 1 . 500 .000 



Fig 1 . Regional map of Iran , location of study area and monitoring stations 
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Fig 2. Annual SPI & RDIst for each hydrological year (Oct-Sep) for yasuj station 




Fig 3. Spatial distribution map of SPI in 2010 
in the study area 



Fig 4. Spatial distribution map of RDI in 2010 
in the study area 



